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Abstract. We consider a model in which a pseudo-scalar field a rolls for some e-folds during 
inflation, sourcing one helicity of a gauge held. These helds are only gravitationally cou¬ 
pled to the inhaton, and therefore produce scalar and tensor primordial perturbations only 
through gravitational interactions. These sourced signals are localized on modes that exit the 
horizon while the roll of a is signihcant. We focus our study on cases in which the model can 
simultaneously produce (i) a large gravitational wave signal, resulting in observable B-modes 
of the CMB polarizations, and (ii) sufficiently small scalar perturbations, so to be in agree¬ 
ment with the current limits from temperature anisotropies. Different choice of parameters 
can instead lead to a localized and visible departure from gaussianity in the scalar sector, 
either at CMB or LSS scales. 
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1 Introduction 

There is currently strong experimental effort to detect the gravitational waves (GW) gen¬ 
erated during inhation. This signal is conventionally parameterized by the tensor-to-scalar 
ratio r, which is the ratio between the GW power spectrum (summed over both polarizations) 
and the density power spectrum at a given large scale. A joint BICEP2/Keck Array and 
Planck analysis has recently reported the 95% CL bound r < 0.12 at the scale 0.05 Mpc“^ 
[1]. It is expected that this limit will be improved by several ongoing and forthcoming CMB 
polarization measurements. For example, a statistical uncertainty a (r) = 0.001 or below is 
quoted in the proposed stage 4 of ground-based CMB experiments [2] . 

It is commonly stated that observing the GW generated during inhation provides both 
(i) the energy scale of inhation, and (ii) a, lower bound on the excursion of the inhaton held 
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during inflation (the so-called Lyth bound [3] , which holds in the case of single field inflation). 
In terms of r, these results read 


yi/4 ^ GeV 



A</) ^ / r \ 1/2 
~Mp ~ VOdil/ 


( 1 . 1 ) 


The observation of an inflationary GW signal does not guarantee that the expressions (1.1) 
are valid, as they rely on the assumption that the observed GW are vacuum modes of the 
metric, amplihed by the inflationary expansion, while it is possible that the observed GW 
are sourced by some other field during inflation. Several recent studies have shown that it is 
rather nontrivial to realize such mechanisms. Whatever the GW source is, it couples at least 
gravitationally not only to the tensor metric perturbations, but also to the scalar ones, and so 
it unavoidably affects also the scalar perturbations to some degree. The source typically leads 
to a statistics that is signihcantly less gaussian than that of the vacuum modes. Therefore, 
respecting the stringent limits imposed by the gaussianity of the scalar perturbations often 
restricts the source term to a too low level to lead to visible GW. 

This problem manifest itself for example ^ in what is probably the most studied mech¬ 
anism of GW generation [24-27]. In this mechanism, the mass of a scalar field x depends on 
the value of the inflaton, and it is arranged in such a way that it vanishes when the inflaton 
reaches a given value 0* during inflation. As the inflaton crosses (/>*, a burst of quanta of 
the field x is generated non-perturbatively [28] and can source GW. In all the implementa¬ 
tions of this mechanism that have been worked out in some details the GW production is 
however negligible; as soon as the inflaton moves past (/>*, the quanta of x rapidly become 
non-relativistic. This suppresses their quadrupole moment, so that the sourced GW signal 
is below the observational level [26] . On the contrary, the production of scalar perturbations 
from these quanta is not suppressed. The non observation of these non-gaussian scalar modes 
puts strong bounds on this mechanism. 

To have a successful mechanism for the production of a visible GW during inflation it is 
therefore crucial to maximize the ratio between the sourced GW and the sourced scalar modes 
[26]. ^ For example, one could imagine that the mass of x is controlled by a field a different 
from the inflaton. The absence of a direct coupling between the inflaton and x decreases 
the amount of scalar perturbations produced by the latter. Even in this case, however, the 
gravitational coupling between x and the inflaton leads to a much greater production of scalar 
than tensor perturbations, due to the non-relativistic nature of x [26]. This led ref. [26] to 
study the case in which the source of GW is both (i) only gravitationally coupled to the 
inflaton, and (ii) a relativistic gauge field A^. 

This is a modification of the mechanism studied in [31], in which the vector field is 
produced by a pseudo-scalar inflaton cj). It was shown in [32] that the pseudo-scalar coupling 

^ Other mechanisms for the generation of GW which are alternative to the standard vacuum production, 
but are still embedded in the inflationary scenario include the use of spectator fields with low sound speed [4- 
6] , the modification of the dispersion relation of the tensor modes in the effective-field-theory approach [7, 8] , 
varying sound speed of the tensor [9], the strong tachyonic growth of chiral tensor modes in the chromo-natural 
inflation [10-12] and in the gauge-fiation model [13], and preheating [14-21] (in this last case the produced 
GW signal is typically at scales much shorter than the GMB ones). It has also been investigated [22, 23] 
whether the presence of many light degrees of freedom could modify the standard relation (1.1) between the 
scale of inflation and r; as shown in [23], this does not appear to be the case. 

^ We note that also in warm inflation [29] the tensor-to-scalar ratio is reduced compared to the standard case 
due to the dissipative effects and the much greater production of scalar with respect to tensor perturbations 
due to the coupled inflaton [30]. 
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jF F ^ amplifies one helicity of the vector field. This helicity mostly sources the GW of 
one definite chirality [31], through the two-body interaction 5A 5A —>• 5g. The chiral 
nature of the signal is particularly interesting, ^ since it leads to violation of parity in the 
CMB that can [35-38] allow to discriminate the GW produced by this mechanism from the 
vacuum ones. ^ However, the vector quanta also inverse decay to inflaton perturbations, ® 
5A + 5A —)• dcj). This again would lead to an exceedingly large GMB non-gaussianity, so that, 
once the CMB limits are respected, the produced GW are undetectable [41, 42]. The scalar 
non-gaussianity can be rednced if the inflation sources a large number N of vector species, so 
that the inflaton perturbations sourced by these vectors are gaussian due to the central limit 
theorem (obtaining a visible GW signal however requires N ~ few hundreds — thousands 
[31]). An observable GW signal may be also achieved if the CMB anisotropies are due to a 
curvaton [31]. Ref. [26] studied instead the case in which the rolling pseudo-scalar is a field 
a decoupled as much as possible from the inflaton (thus realizing the conditions (i) and (ii) 
mentioned above ^). It was shown there that, in this case, the 5A -|- 6A —)• 64> gravitational 
interaction produces a negligible amount of inflaton perturbations, so that the mechanism 
can indeed lead to a visible GW signal. 

The mechanism of [26] was further studied in [45], as a possible way to generate a blue 
GW spectrum, as required to reconcile the original inflationary BICEP2 interpretation [46] 
of their signal, and the upper bounds on r at larger scales. Some specific potentials for a 
were considered in [26], and it was noted there that for the mechanism to be successful, one 
should also ensure that a does not contribute to the observed CMB perturbations. It was 
mentioned that the simplest way to achieve this is to impose that a becomes massive short 
after the CMB modes leave the horizon. This issue was further studied in [47] and in [48], 
where it was shown that the 6(j) production is dominated by their linear interaction with the 
6a quanta that are sourced by the gauge field and that the amplitude of the 64> fluctuations 
sourced this way is proportional to the number of e-folds during which a is rolling. In this 
work we note that the 64> fluctuations are sourced only for modes that leave the horizon when 
a ^ 0 and that the constraints from non-observation of a nonvanishing scalar bispectrum are 
relatively weak for modes with I < 100 or so. As a consequence, if d ~ 0 when modes with 
I ^ O(IO^) leave the horizon, then the constraints from non-gaussianity can be more easily 
evaded. 

In the computations of [26] and [47], a was assumed to be constant. Here we study the 
generation of primordial perturbations in a more realistic setup, namely in a model in which 


^In this expression, / is a mass scale, often denoted as the “axion scale”, is the field strength of 

the vector, and F^'' = Y^^=Fap its dual. The quantity is totally anti-symmetric, and normalized to 

^0123 _ 

^Moreover, the (f>FF interaction does not lead to such a ghost instability as the mechanism of chiral GW 
production from a gravitational Chern-Simons term coupled to the inflaton [33, 34]. 

®The phenomenological constraints on the production of chiral magnetic fields and their subsequent sourc¬ 
ing of chiral gravitational waves have been studied in [39, 40]. 

®At the technical level, in the presence of a direct coupling the sourced scalar power is enhanced with 
respect to the sourced tensor power by a factor of cf. eqs. (8) and (10) of [41] and the discussion in section 

4.1 of [31]. 

similar mechanism that also satisfies these conditions, and that also leads to observable GW, has been 
recently studied in [43]. In this work, the sourcing vector field is produced by a rolling dilation field cr (only 
gravitationally coupled to the inflaton) through the Ratra / (a) F^ mechanism [44] . 
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the field a evolves in a given potential. The simplest potential for an axion held, 



( 1 . 2 ) 


allows for a very small a both at early and late times, when a is, respectively, close to the 
maximum (a = 0) and the minimum of the potential (a = /vr). The speed has a peak 
at an intermediate time, when a is between the maximum and the minimum (we assume 
that Va is subdominant with respect to the inhaton potential; moreover, we assume that the 
inhaton potential is very hat, so that the Hubble rate H can be treated as constant; this is 
the most interesting regime for our discussion, since a very hat inhaton potential corresponds 
to unobservable vacuum GW). The peak lasts for a number of e-folds roughly of O 
where m is the curvature of (1.2). So, remarkably, the simplest axion potential is a perfect 
candidate for generating a visible GW signal, while keeping the 6(f) production under control. 
In this work we show that this is indeed the case through explicit computations. 

We present some specihc examples (namely, some choice of parameters in the model) 
for which the sourced tensor mode strongly dominates over the vacuum one at large scales, 
leading to observable B modes of the GMB polarization. The GW signal also leads to a 
marginally observable TB correlation (as a consequence of the broken parity invariance of 
the mechanism) and to a well observable (high signal-to-noise ratio) BBB correlation. At 
the same time, in such examples we find no statistically significant signatures in the TT, and 
TTT temperature correlators. 

The plan of the work is the following. In Section 2 we present the model, the background 
evolution, and the vector field production. In Section 3 we study the cosmological pertur¬ 
bations (scalar and tensor modes) sourced by the vector field. In Section 4 we summarize 
our results for the two- and three-point scalar and tensor perturbations, and we discuss their 
phenomenology. In Section 5 we present our conclusions. The work is complemented by six 
appendices. In Appendix A we compute (in WKB approximation) the gauge field produced 
in the case of nonconstant a. In Appendix B we review the computation of scalar modes 
produced in the case of constant a. In Appendices C and D we give details of, respectively, 
the scalar and tensor mode computation. In Appendix E we present some properties of the 
bispectra produced in the model. In Appendix F we estimate the departure from gaussianity 
of the statistics of the sourced modes. 

2 Model, background evolution, and vector field production 

We will consider a system containing the inflaton cf) and a second rolling field a which interacts 
with the 17(1) gauge field via an axionic coupling, so that the lagrangian reads 


£ = _1 ( 5^)2 _ 1 ( a ^)2 _ _ 1^2 _ 


( 2 . 1 ) 


The rolling of a provides a time dependent background for the gauge field and amplifies 
its vacuum fluctuations. Such a phenomenon, on a de Sitter Universe with expansion rate 
H, is controlled by the dimensionless quantity 


^ 2fH ’ 


a a 


which must be larger than unity or so to give a significant effect [32]. 
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In order to decouple as much as possible the fluctuations in the gauge field from the 
scalar perturbations that are strongly constrained by CMB measurements, we will assume 
as in [26] that there is no direct coupling between the inflaton and the field a, and we will 
write a) = V^{4>) + I4(fT). 

As we have discussed in the Introduction, if a is approximately constant throughout 
inflation, then an exceedingly large non-gaussian component of the metric perturbations is 
generated. However, a large weight in the constraints on the primordial bispectrum from 
cosmic microwave background radiation measurements is carried by relatively high multi¬ 
poles, t = 0(100 — 1000), whereas a possible observation of primordial tensor modes relies 
on measurements at larger scales, corresponding to ^ ~ 0(10 — 100). As a consequence, we 
will focus on the case where a is significantly different from zero only during the epoch when 
scales of the order of 10“^ — 10“^ times the current size of the observable Universe have left 
the inflationary horizon. 

Remarkably, the most natural potential for an axion-like field such as u, namely eq. 
(1.2), does the job: assuming <C U^, one finds the slow roll solution 


a = 2f arctan 




(2.3) 


where 5 = A'^/(6ff^ /^) and where we have denoted by the time at which a 
consequence. 


fH5 


cosh [5H (t — U)] 


The slow roll condition gives 


I /. As a 
(2.4) 


a 

3H& 


5 

3 


tanh [6H {t — U)] , 


(2.5) 


so, to ensure that the slow-roll solution is valid, we require that <5 <C 3. Finally, the parameter 
^ is given by 


a& _ 

2Hf cosh [(5H (t — U)] 



( 2 . 6 ) 


where and by a* are, respectively, the value of ^ and of a at t = t*; we note that = a 5/2. 
The value of ^ is significantly different from zero and of the order of only for (2^*)“^/'^ < 
a/a* <(2e,)V'5. 


2.1 Amplification of the fluctuations of the gauge held 

The last term in the lagrangian (2.1) is responsible for the amplification of the vacuum 
fluctuations of the gauge field. To see this, we choose the Coulomb gauge and we define 


Ai{T, x) = 


(fk 




(27r 


, 3/2 


E 

A=± ' 


d^k 


(2^: 


, 3/2 




(fc) Ax{t, k)dx(k 


e*^'^ + h.c. 


(2.7) 


where the helicity vectors satisfy k ■ =0, ik x = ±ke^^\ = 1 and 

e(±) . e(±) = 0. Explicitly, if = (sinScosiji, sin0sin(/), cos 6), then 


(^) ~ ^ (/>, cos 6 sin 0 ± / cos 4>, — sin 6) . (2.8) 
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The mode functions A± satisfy the equation 


Al+ [k^ =fk 


a a 


f 


A± = 0, 


(2.9) 


where the prime denotes differentiation with respect to the conformal time r. In the case of 
spatially flat, inflating Universe with Hubble parameter H and scale factor a(r) = —1/{H r),® 
the second term in parentheses in eq. (2.9) reads kaa'/f = —2k^/T. 

For constant values of eq. (2.9) can be solved exactly [32]. In our case, we must solve 
the following equation 


A'i + k^± 


4k^^ 


(r/r^f + (n/rf 


A± = 0, 


( 2 . 10 ) 


where is defined in (2.6), and r* = —l/(a*H). Without loss of generality, we set > 0. 
As a consequence, only positive helicity modes are amplified, and we can neglect the A_ 
mode. 

The solution of eq. (2.10) cannot be written in closed form. However we can obtain a 
good analytical approximation of the time dependence of A_|_(r, k) using the WKB approxi¬ 
mation, finding (see Appendix A) 


A+ (r, k) ~ 


—r 


8fc^(r) 


1/4 


where we have defined 

A (r, k) = N X*, d] exp 


A (r, k) , (r, k) ~ 

1/2 


4^: 


1 + (5 V F 


-2t 

( 5/2 


-I 1/4 


A{T,k) , 


— i-kr) 


1/2 


( 2 . 11 ) 


( 2 . 12 ) 


with X* = — A:r* = k/{a^:H). We determine the time-independent normalization factor 
N X*, (5] by matching A_|_ at late time —kr <C 1 to the full numerical solution of (2.10). 
We choose the arbitrary initial phase of A+ such that N [(^*, x*, d] is real and positive. As a 
result, in what follows we will use the following expression for the gauge field 

,1..^ /^\ a+ ( —k ) , (2.13) 


Ao = 0, Ai[T,k] = 


{2tt 


_ ik-x (+) 

, 3/2 i 


with the function A+(r, k) given by eq. (2.11). 

For future reference we give here also the “electric” and “magnetic” field, ® which are 
related to the vector potential by 


E- — _ L 4' 

■CJi — c) 1 


1 


Bi — 2 U/fc djAf^ , 


(2.14) 


and therefore read 

Ei (t, k^ = - J 

Bi (T,k] = 


Sk 

(27r)^/^ 

d?k 


^ik-x ^(-1-) 


(.) 




fc^(r) 


-11/4 


(27r 


, 3/2 




-2t 

_fc3.ll/4 


L8^(r)J 


A (r, k) a+ yij + y—k 
A (r, k) a+ (k^ + at (—k 


(2.15) 


^In this work, we disregard all terms that are subdominant in a slow-roll expansion. 

®We do not need to identify the gauge field considered here with the Standard Model photon. Nonetheless, 
we sometimes use electromagnetic notation for convenience. 
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using the approximate solution (2.11). 


3 Sourced primordial perturbations 

This section is divided in two subsections. In Subsection 3.1 we study the scalar modes 
generated by the vector modes computed in the previous section. In Subsection 3.2 we study 
the sourced tensor modes. 


3.1 Generation of scalar perturbations 

The action (2.1) contains two scalar and two tensor degrees of freedom as dynamical modes. 
After taking the spatially flat gauge, the scalar sector metric reads 

(r) [— (1 + 2cj)) dr'^ + 2diB dx^dr + dijdx^dx^] , (3.1) 

and, after solving for the non-dynamical variables (j) and B, we decompose the remaining 
physical modes (/>(r, x) and cj(r, x) as 


<^(r,f) = 4){t) + 


(T(r, f) = fj (r) + 


d^k 


(27r 


,3/2 


d?k 


Ak-x 


Ak-x 


a(r) 

Qa (k 


(27r 


,3/2 


(3.2) 


and denote { 4 >i,(j) 2 ) = (</>,o') and (^Qi,Q 2 ^ = 
then reads 

4ee [Q^] = 11 drdA [g't Q/ _ gt 

where prime denotes derivative with respect to the conformal time r, and 


a(r) 

The free part of the action for Qi 


— - 6ij + a Vii + ( 3 — 




A A A ^'A'i , 


^ \ Z2Z1 , (A'y. , A/v ) 

2M| a'2 ) M2 + M^a' 


with I/j = dVld4>i. The interaction term in (2.1) gives 


a 


a 


Sint = - j = j dAaAjE-B , 


(3.3) 


(3.4) 


(3.5) 


where the “electric” and “magnetic” fields are defined in (2.14). The equations of motion for 
Q(/) and Q^, derived from (3.3) and (3.5), are 


+ M^^Qa- — 0 , 

(1^+*^' + *-) <3. + J 


(3.6) 


In principle the inflaton also conples to the gauge sector through gravity and receives contributions of the 
type 5A + 5A -A 5(j>. This coupling term is a Planck suppressed operator, and parametrically one may expect 
its effects to be as large as the ones from 5^4 + 5A -a 5a ^ S(j>, which is of our main interest in this work. 
However, 5a is produced near horizon exit and keeps sourcing 5<j) while being outside the horizon, whereas 
the direct production of from the gauge field occurs only near horizon crossing and, as shown in [26], is 
negligible. Therefore the latter mechanism, mediated by 5a, provides the leading contribntion from the gauge 
field to the curvature perturbations. 
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Expanding Mjj to first order in the slow roll parameters, we have 




2 + 9 + 3 


3 

2 + 9 €(7 + 3 



where we have defined the slow-roll parameters as 


(3.7) 


^<i>i 


2M^m ’ 


^4>i — Y ’ 


(3.8) 


with no summation for the repeated i indices. In the following computation, we focus on 
the production of da from the gauge field 6A and its subsequent sourcing of 5(j), namely the 
process 6A + 5A 5a ^ 54>. In this regard, we only take the dominant terms in (3.7) and 
approximate (3.6) as 






— 2 's/Qo 

(Pi 


d"^ 2 2 \ - a? . — 


Q-ik-x B = 


= 5,7 (r, fe) 


(3.9) 

(3.10) 


We note here that we are disregarding the backreaction of the quanta sourced by on 
the sourcing Qo- modes; this effect is of higher order in slow-roll. Let us also emphasize that 
since we are interested in the feature due to the change in a during the trajectory, we do not 
neglect the time dependence of e^, while we assume constant e^p. 

In the spatially flat gauge, the scalar curvature perturbations are related to the inflaton 
perturbations by 


C (r. fc) =. -1 i0 (r,«) = Q, (r. t) 


(3.11) 


assuming <j) > 0. To find Qp) ^r, kj , we solve eq. (3.9) by separating Q</) into two parts, 

= + , (3.12) 

where is the homogeneous solution of eq. (3.9) and is its particular solution. These 

two contributions are statistically uncorrelated. We decompose the operator corre¬ 

sponding to the vacuum fluctuations, as 

= Q^^\T,k)a{k) -hQ®*(r,A:)al( - fe) , (3.13) 

a general two scalar field model, the curvature ^ is a linear combination (p (t) = ci (t) S4> {t) + C 2 (t) 3a (t), 
where the coefficients ci ,2 depend on the background (the orthogonal combination being an isocurvature mode). 
With the choice of potential (1.2), and with our assumptions on the parameters, the field a becomes massive 
short after the CMB modes are produced. As cr becomes a massive field in an inflationary universe, its energy 
density and pressure very rapidly drop to zero, and so does C 2 . The only potentially observable effect of 6a 
is through its couplings to the inflaton and metric perturbations, which are effective only as long as ct yf 0. 
The dominant among these interactions is the linear coupling to Scj) that we are computing here (all the other 
couplings are nonlinear in Sa, and highly subdominant; they can be disregarded, once we impose that the 
effects of the linear coupling are below the observational level). 
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where and a are the creation and annihilation operators for respectively, and the 

mode function Q® is solved to be, with the Bunch-Davies initial conditions, 




—ikr 




1 - 


kr 


(3.14) 


The particular solution ^ is obtained by solving eq. (3.10) and plugging its solution 
into (3.9): 


= J dr' Gk{T,T') ^ j dr” Gk{T',T")Sa{T",k) , (3.15) 

where the retarded Green function Gk{T,T') reads 

Gk{T,T') = 0(r-r') [j3/2(-fcr)y3/2(-/cr') -T3/2(-fcr) J3/2(-/cr')] , (3.16) 

where J and Y denote the Bessel function of real argument. In this work for the first time we 
evaluate the sourced solution (3.15) in the case in which a is time-dependent. In Appendix 
B we show that the solution reproduces that of [47] (up to the correction of a typo) in the 
case of constant a. 


3.2 Generation of tensor perturbations 

Let us next focus on the tensor perturbations, considering a metric of the form 


ds'^ = (r) —dr'^ + ( 5ij + hij (r, x) j dx^dx^ 


(3.17) 


where is transverse and traceless. The leading expression for hij is obtained by expanding 
the action to second order in hij, including the first order interaction term with the gauge 
field: 








^ij dij 





(^Ei Ej -\- Bi Bj^ 


(3.18) 


To obtain a canonically normalized field describing tensor modes in Fourier space we 
decompose 



2 

Mp a{T) 


d^k 


Ak-x 


(2vr) 


3/2 


E 




k] Qx (r. 


A—^ 



(3.19) 


where the polarization operators are 

± (k) = (k) ef (k) . (3.20) 

We neglect the scale dependence of the homogeneous solutions of eq. (3.6) that can be induced by relatively 
large values of e„ and rj„-, we have checked that the mixing angle, see eq. (B.2), between and can be 
made small enough that such a scale dependence will affect only the unobservable Qa mode and will not leak 
into the metric perturbations, that are associated to Q 4 ,. 

^^Interaction terms of the form hhFF give contributions to the two-point function of the tensor perturba¬ 
tions that are of the same order in as those of the form hFF. However the terms of the form hFF give 
a contribution to the tensor power spectrum that is proportional to , much larger than those of the form 
hhFF, whose contribution is proportional to which therefore will be neglected [49]. 
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The equations of motion for Qx 






(27r)3/2 


^ ^-ik-x 


Ei Ej + Bi Bj 


= Sx{t, k ), 
(3.21) 


are solved, as in the scalar case considered in the previous subsection, by separating Qx into 
a vacuum mode solution of the homogeneous equation, and a sourced mode (also 
in this case the two modes are statistically uncorrelated). The vacuum mode is given by 


QT (^) = (r, k) OA (fe) + h*x (r, k) dl [-k] , 

^—ikr 




hx (r, k) = 




k T 


(3.22) 


where creates gravitons of helicity 2 A. 

For the sourced mode, we have the formal solution 

= j dr' Gk (r, r) 5a(t, k) , (3.23) 

where the retarded propagator is the same as that for the scalar perturbations, eq. (3.16). 


4 Results and phenomenology 

The gauge quanta produced by the rolling pseudo-scalar a source quanta of a through inverse 
decays 5A + 5A —)• 5a. They also source inflaton perturbations and gravity waves through 
the 2 —)• 1 processes 5A + 5A ^ 64> and 6A + 5A —)• hx- Under the assumption of no direct 
coupling between the inflaton and a in the potential, the 5A + 5A —?• Scj) interaction is of 
gravitational strength, and negligible with respect to the gravity wave production [26]. For 
this reason we disregard it in this paper. On the other hand, the inverse decay into quanta of 
a can produce a large signal, which has a much stronger departure from gaussianity than the 
vacuum mode [41]. In the model under consideration the field a has a completely negligible 
energy density at the end of inflation / beginning of reheating (both at the background 
and at the perturbation level), so that the inflation perturbations can be identified with the 
observed curvature perturbations C) see eq. (3.11). However, 5a is not completely irrelevant 
for observations: as long as a is rolling (see eq. (3.7)), there is a linear gravitational coupling 
between perturbations of a and of the inflaton, leading to conversion of some 5a into 54> 
quanta [47]. 

In the previous Section (see also appendices C and D) we computed the curvature 
perturbations ([ and the gravity waves produced in the model. The results are summarized in 
Subsection 4.1. In the following Subsection 4.2 we study the phenomenological implications 
of these results. 

4.1 Scalar and tensor correlators 

We are interested in the power spectrum and bispectrum of the scalar curvature C and of the 
gravity wave polarizations h±. They are defined as 

(c(fc) ({k'))^^Aj,^(k)sm(k + k') , 

(C(tl) C(fe) ([k3)')sB((ki,k2,k3}S'^>(kj + k2 + h) , (4.1) 
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and analogously for h±. The scalar curvature C and gravity waves h± produced in the model 
(2.1) are the sum of a vacuum component and a component sourced by the gauge quanta. 
The two components, which we denote, respectively, as 


C = + , h± = hf + h^l\ 


(4.2) 


are uncorrelated among each other, and therefore the power spectra and bispectra of these 
modes are the sum of the vacuum and of the sourced spectra: 


V^{k) = vf{k) + v^l\k) , B^{k) = Bf\k) + B^l\k) , 

and analogously for h±. 

The vacuum power spectra are given by the standard relations [50] 


Vf {k) 




(k) = (k) 


7r2M2 ’ 


(4.3) 


(4.4) 


where we disregard subleading corrections in slow-roll. This gives the vacuum relation 


^vac 


'pW _j_ 'p(^) 


r, 


( 0 ) 


~ 16e^ 


(4.5) 


for the tensor-to-scalar ratio. The vacuum bispectra are known to be well below the current 
observational limits, and we disregard them in this work. 

The sourced terms and present a localized bump in momentum space. This 
corresponds to modes that leave the horizon when & is close to its maximum value. The 
gauge field production is exponentially proportional to the parameter given in eq. (2.6) 
[32]. Therefore the height of the bump is exponentially proportional to the maximum value 
acquired by The width of the peak decreases with increasing values of the parameter S. 
This is because, as we have see in Section 2, the field a has a significant evolution only for a 
number of e-folds AN ~ j, and, therefore, the larger 6 is the fewer are the modes produced 
while a has a substantial rolling. 

As we show in Appendices C and D, we can factorize the following parametric and 
functional dependence from the scalar and tensor spectra 


) (fc) 

r? {k) 

B( {ki, k2, ks) 

Bx {ki, k2, ks) 




(■(p'P^ ^k) /2,A ( 1 ) ~ ) 


1 2 




1 3 


l2 u2 u2 
^2 ^3 


[ T ) I ) I ) 0 ) ) 

Ay* A/j- A/* 


(k) 


1 3 


t.2 u2 u2 

fi,l ^3 


, ki k2 ks 
•'SA [ s*) 

Ay* Ay* Ay* 


A — , 


(4.6) 


For brevity, we denote the functions at the right hand side as fij, where i G {2, 3} and j G 
{(C, +, —}. These functions are dimensionless. The functions encode the full dependence 
of the bispectrum on the momenta ki. In eq. (E.2) we provide an approximate relation 
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Figure 1. Comparison of the exact function / 2 ,^ (green dots) and of the approximate form (red 
dashed line) given in eq. (4.7), for the two values of 6 that we have studied and for one illustrative 
choice of The function f 2 Xi defined in (4.6), depicts a bump in the sourced scalar power spectrum, 
featuring a short period of a relatively fast rolling of a and the production of the gauge field during 
this stage. A bump feature is also present in the tensor sector. The phenomenological consequences 
of these localized signal are studied in Subsection 4.2. 


for the shape dependence, written in terms only of the two point function and of the three 
point function in the exact equilateral case. The expression (E.2) is exact by construction 
on equilateral triangles, and it is very accurate where the signal is maximum (see Figure 
10). Therefore, to study the phenomenology of the scalar and tensor two- and three-point 
correlation function we simply need to provide the functions, f 2 j, as well as /sj for equal 
momenta. 

We first studied the momentum dependence of these functions for hxed values of and 
6. We found that they are well approximated by 


fi,j ( , ) 


fi,j exp 


In' 




(4.7) 


Namely, the momentum dependence is encoded by the three functions a, which in 
turn depend on the evolution of & through the model parameters and 6. The function (4.7) 
has a bump at k = where it evaluates to The parameter controls the width of 

the bump. For each choice of and 5, the values of /T, x^ and aij are obtained by fitting 
the right hand side of (4.7) to reproduce the position, height, and curvature of the bump. In 
Figure 1 we verify the accuracy of the approximate relation (4.7) by showing the comparison 
between the exact and approximate form of f 2 x for two choices of parameters (the chosen 
values have no particular relevance, and have been chosen for illustrative purpose only). An 
analogous level of accuracy is obtained for the other fij functions and for the other choices 
of that we have made. 

Secondly, we need to specify the dependence of /A, on the model param¬ 

eters and 6. We evaluated these functions for 5 = 0.2 and 6 = 0.5, and for = 
{2.5, 3, 3.5, 4, 4.5, 5, 5.5, 6, 7}. We found that the dependence on is rather smooth, and 
can be well described by a second degree polynomial. We provide the fitting functions in 
Tables 1 and 2. As we see from the examples shown in Figure 2, the fitting functions are 
very accurate. 

^^For the three point functions (i = 3), the approximation (4.7) refers to the exact equilateral case, ki = 
k2 = ks = k. 
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{h i} 

In \fu\ ^ 

<7 - 

(Ti,j ~ 

{ 2 , a 

-5.60 + 10.U* -b 0.0947^^^ 

3.26 + 0.435^* +0.0109 

1.41 -0.166^* +0.00962^;^ 

{2, +} 

-7.98 + 10.0 +0.0979 

5.45 + 0.455^* +0.0316 

1.38 - 0.178 C* + 0.0103 C* 

{2,-} 

-13.8 + 9.96 +0.104^;^ 

2.39 + 0.129 C* +0.0214 

1.44-0.169^* + 0.0102 

{3, a 

-4.72 + 15.U* +0.142^;^ 

3.37 + 0.353 C* +0.0142 

1.13-0.143 C* +0.00830 

{3, +} 

-7.77+ 15.U* +0.147^;^ 

5.24+ 0.36U* +0.0360^;^ 

1.11 -0.150 C* +0.00880^;^ 


Table 1. dependence of the functions entering in (4.7), for S = 0.2. Among the entries in the first 
column, only /| _|_ is negative, while the other are positive. 


{h j} 

In 1/^,1 ^ 

- 

(TiJ ^ 

{2, a 

-6.47 + 9.04^* +0.05864^ 

1.64 + 0.630 C* +0.00738^)^ 

0.823 - 0.0872 C* + 0.00558 C'i 

{2,+} 

-6.85 + 9.05^* +0.05964^ 

2.70 + 0.896 C* +0.0187^^ 

0.768 - 0.0993 C* + 0.00608 

{2,-} 

-12.5 + 8.97 C* +0.0656^*^ 

1.22 + 0.396 C* +0.00976^*^ 

0.858 - 0.0813 C* + 0.00530 

{3, a 

-5.98 + 13.6 C* +0.086U*^ 

1.71 + 0.569^* +0.00972^;^ 

0.641 - 0.0792 C* + 0.00495 

{3, +} 

-6.03 + 13.6 C* +0.0870^*^ 

2.63 + 0.841^* +0.0193^^ 

0.606 - 0.0844 C* + 0.00520^)^ 


Table 2. dependence of the functions entering in (4.7), for 6 = 0.5. Among the entries in the first 
column, only /| _|_ is negative, while the other are positive. 




Figure 2. dependence of the fitting / 2 ,^, and cr2,c entering in (4.7), for S = 0.2 (an analogous 
agreement is obtained in the S = 0.5 case, and for th other {i, j} functions). The red dots denotes 
the values obtained by evaluating the power spectrum at = {3, 3.5, 4, 4.5, 5, 5.5, 6, 6.5, 7}. The 
solid lines are the polynomial fits reported in the first row of Table 1. 


4.2 Phenomenology 


In this subsection we study the phenomenology of the model (2.1). The total scalar and tensor 
modes produced during inflation are the sum of a nearly scale invariant vacuum mode, plus a 
peaked sourced signal (see Figure 1 for the scale dependence of the latter). The sourced scalar 
and tensor mode modify the tensor-to-scalar ratio from its vacuum value (4.5), ryac — 16e,^, 
to 


., vf (k) + (k) + yX' (k) 1 + a rf (fc) h+ jk) 

where we have disregarded the sourced mode (given that f^- ^ / 2 +)- Both the numer¬ 
ator and the denominator of the fraction appearing in the last expression of eq. (4.8) have 
been written as 1 + Psourced/^vacuum, so that we can immediately compare the impact on r 
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Figure 3. Black-solid lines: tensor-to-scalar ratio (4.8). Orange dashed (respectively, blue dotted) 
lines: ratio between the sourced and the vacuum tensor (respectively, scalar) power spectrum. All 
ratios are evaluated at the peak of the sourced GW power spectrum. 


of the sourced tensor modes vs. the sourced scalar modes. We see that, relatively to the 
sourced scalars, the sourced tensors are more relevant at smaller values of e^. This is the 
regime of greatest interest for our study, since it corresponds to a small rvac- 

This is confirmed by the contour lines shown in Figure 3, where we show the tensor-to- 
scalar ratio (black solid lines), the ratio between the sourced and the vacuum tensor power 
spectrum (orange dashed lines), and the ratio between the sourced and the vacuum scalar 
power spectrum (blue dotted lines), evaluated at the peak of the sourced GW signal. The 
ratios are shown as functions of and e^. The sourced modes are negligible at the smallest 
values of shown in the figure, so that r ~ rvac — ISe^ there (so that the contour lines of 
equal r are independent - and, therefore, horizontal, in this region). At greater values of 
and, particularly, at the smallest e<^ shown, we see that r rvac- 

Figure 3, and several of the following figures, are comprised of two panels. In the left 
(right) panel we show results for 5 = 0.2 {6 = 0.5). The roll of a is substantial for a number 
of e-folds AN ~ cf. eq. (2.4). Increasing 6 therefore decreases the amount of time during 
which a is rolling. This decreases the amplitude of the produced signal (encoded in the 
parameter cf. eq. (4.7)), as well as its width (encoded in aij). The latter effect is due 
to the fact that, the smaller AN is, the fewer are the modes that exited the horizon while a 
was non-negligible. By comparing the results given in Tables 1 and 2, we can also see that 
increasing 6 decreases more the scalar production than the tensor one. This is also visible 
by comparing the left and the right panel of Figure 3, as well as from the following figures, 
where we can notice that greater values of d corresponds to a greater production of tensor vs. 
scalar modes. This occurs because the inflaton perturbations are sourced by the 6cr modes 
only as long as (t 7 ^ 0. So, a decrease of AN affects more the sourced scalar modes than the 
tensor modes because it decreases both (i) the number of modes that are sourced (namely, 
the width of the peak), and (ii) the interval of time during which the 6cr modes (produced 
by the gauge field) can be converted into inflaton perturbations. 

In the following analysis, we show that the model admits choices of parameters that 
result in a sourced signal visible in correlators involving the B-mode of the CMB polarization, 
and that are consistent with the well-analyzed TT and TTT correlators. To have visible 
effects from the sourced GW, we choose values of fe* leading to a peak of the sourced signal 
at large GMB scales. Specifically, we choose ranging from 7 x 10“® Mpc“^ (affecting only 
the first few multipoles) to 5 x 10“^ Mpc“^ (affecting multipoles up to the first acoustic peak). 
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Figure 4. Solid lines: Largest value of (controlling the amount of produced quanta) allowed by 
the WMAP TT data, obtained as explained in the text, as a function of e^, and for different values of 
fc* (controlling the scale of the bump of the sourced modes, see the next figure). Dashed lines: ratio 
of ^at the peak of the GW bump. 

4.2.1 CMB power spectra 

The main goal of this work is to provide a concrete example of a particle physics process that 
can enhance the GW signal from inflation, while being “consistent with the TT and TTT 
data”. With this, we mean that this additional signal should not significantly worsen the fit 
to the TT and the TTT data of the standard cosmological model which does not include this 
signal. The signal that we are studying manifests itself at the largest cosmological scales; we 
compare it against the latest WMAP data [51, 52], which are cosmic variance limited at such 
scales. We fix all the cosmological parameters consistently with the best fit values reported 
in [52], denoting this set of values as Co; we denote by To the likelihood of the fit to the 
WMAP TT data of the standard cosmological model with these values of the parameters: 

Co = . (4.9) 

We then add the sourced scalar and tensor modes obtained from the mechanism presented in 
the previous sections. For illustrative purposes, we fix three values of fe* (controlling the scale 
at which the sourced signals are peaked) and two values of 5 (controlling the width of the 
sourced peak). For any of these six choices, we then increase (the parameter that controls 
the amplitude of the sourced signal) until the likelihood of the fit decreases by a factor e^ 
with respect to Co, namely, we find the value = G nmu . for which + 

experimentalist fitting the WMAP data with the cosmological parameters Cq, would obtain 
a fit that is two sigmas worse if he/she includes the sourced signal instead of the standard 
cosmological model. In fact, the value obtained in this way is conservative, as we 

do not vary any of the cosmological parameters Cq. Varying them, we can only improve 
the fit of the WMAP data for that given value of and so likely obtain larger values for 
limit- However, changing values in Co may give a disagreement with data at smaller angular 

Specifically, we fix the physical baryon density = 0.02264, the physical cold dark matter density 

Dc = 0.1138, the ionization optical depth r = 0.089, the Hubble constant Ho = 70.0 km/sec/Mpc. The 
parameters relevant to the initial condition, i.e., the scalar amplitude = 2.41 x 10“®, the scalar spectral 
index Us = 0.972 and the WMAP pivot scale ko = 0.002 Mpc“^, fix our vacuum power spectrum = 

A^{k/ko)"‘'~^■ We have also assumed a flat universe, with 3.046 relativistic species and no massive neutrino. 
^®We use the WMAP power spectrum likelihood code: http://lambda.gsfc.nasa.gov/. 
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Figure 5. Temperature-temperature CMB power spectrum. The final WMAP data [51, 52] are 
compared against the theoretical curves, evaluated for = 10“® and for 6 = 0.2 (left panel) and 
5 = 0.5 (right panel). In each panel we show the theoretical curves for three different values of fc* 
(corresponding to three different locations of the peak of the sourced signal) and for the limiting value 
i* = limit, obtained as explained in the text. 

scales than the WMAP ones, and for this reason we do not vary such parameters. As we 
show below, in most cases the values of ^*,iimit obtained with this procedure already provide 
a visible GW signal, which is the goal of this present analysis (in summary, we are not 
interested in providing precise Bayesian limits on within this model - for which we should 
provide priors, marginalize over all the other parameters and include smaller scales data - 
but only in the goal specified at the beginning of this subsection). 

In Figure 4 we present the value of obtained with this procedure as a function 

of the slow roll parameter e^. In the left (right) panel of the figure we set S = 0.2 (0.5), 
corresponding to a significant a for about 5 (2) e-folds. In each panel we fix fe* = 7 x 
10“® Mpc“^ , /c* = 5x 10“^ Mpc“^ , and = 5x 10“^ Mpc“^, producing a bump, respectively 
at the very largest angular scales {£ ^ 5), on the rise of the hrst peak, and in the region 
around the first two peaks. In the same figure we also show with black dashed lines the 
ratio at the peak of the GW bump. We see that the allowed amount of scalar 

signal strongly depends on the scale. In the examples with the bump at the largest scales, 
the sourced signal can be as large as the vacuum one at the peak, due to the large cosmic 
variance present at those scales. A significantly smaller fraction, O (1% — 10%), is allowed in 
the examples in which the signal affects the acoustic peaks. 

By comparing the left and right panel of Figure 4 we again see that, for any fixed value 
of the sourced signal is stronger at small values of 5. 

In Figure 5 we show the TT power spectrum obtained for the same choices of 5 and fe* 
as in Figure 4, for = 10“®, and for the corresponding value of limit- The theoretical 

curves present a bump due to the sourced scalar modes. As we already mentioned, the bump 
ranges from the lowest i multipoles (for the smallest /c* chosen) to ~ the first two acoustic 
peaks (for the largest chosen). 

In Figure 6 we show the self-correlations between the B modes of the GMB polarizations 
(BB) sourced by the tensor modes when limit- The parameters e^, 6, and k^ are chosen 

as in the two previous figures. For comparison we also show with black dashed lines the BB 
correlation obtained for a scale invariant r of 0 . 1 , 10 “^, and 10 “^ (from top to bottom, 

verified that, apart from the {<5 = 0.5, fc* = 7 x 10“® Mpc“^} case, the sourced GW give a negligible 
contribution to the TT signal, and in all the other cases the limits shown in the figure are due to Csourced- 
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Figure 6. First row: Forecasted signal-to-noise ratio for the detection of the B-mode auto power 
spectrum in a realistic CMB experiment with Planck-\\kB sensitivity (Colored solid lines), and in a 
cosmic variance limited (namely ideal noise-free) CMB experiment (Black dotted line), as a function 
of the maximum multipole t included in the analysis. Second row: coefficients. The model 

parameters are chosen as in the previous figure. 


respectively). As we already mentioned in the Introduction, the proposed stage 4 CMB 
experiments claim an expected statistical uncertainty a (r) = 10“^ or below [2] in the scale 
invariant case. If this is achieved, the theoretical curves chosen in Figure 6 appear to be within 
observational reach. We recall that = 10“® corresponds to a vacuum tensor-to-scalar ratio 
I'vnc — 1-6 • 10“^. Therefore, the enhancement of the BB signal visible in the figure is entirely 
due to the sourced tensor modes. The enhancement is present at progressively larger i for 
increasing values of /c* shown (namely for bumps of gauge field production at progressively 
smaller scales). By comparing the left and the right panel of Figure 6 we observe that BB can 
reach greater values at increasing 6. This is consistent with what we have already mentioned: 
at fixed both the sourced scalar and tensor modes decrease with increasing 5. However, 
the scalar mode decreases more. Therefore, at larger values of 6, larger values of can be 
compatible with the WMAP TT bounds (cf. Figure 4). Such values lead to a larger amount 
of sourced tensor modes. 

The signal-to-noise (S/N) ratio shown in the figure is evaluated through 

(sy 21+1 (ci“>\ 

Here, (C^^)^ corresponds to the signal given by our theory, while the other terms in this 
relation account for the uncertainty of the BB power spectrum in a given experiment. For 


(4.10) 
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Figure 7. First row: Forecasted signal-to-noise ratio for the detection of the TB correlation in a 
realistic CMB experiment with Plancfc-like sensitivity (Colored solid lines), and in a cosmic variance 
limited, CMB experiment (Colored dotted lines), as a function of the maximum multipole t included 
in the analysis. Second row: TB spectra. The parameters are chosen as in the previous figures. 


simplicity, we here (and also in the other S/N estimations) assume a full-sky isotropic CMB 
measurement, thus, the summations in terms of m disappear in the S/N formula. The data 
spectrum in a given experiment is regarded as the sum of the signal and instrumental noise, 
reading = Cf^+N^^. In the paper, we analyze two different types of measurements: a 

realistic measurement including a P/ancA:-level noise spectrum [55] (as described in Appendix 
A of [56]) and an ideal noise-free cosmic variance dominated measurement (i.e., = 0). 

The results are shown in the upper panels of Figure 6. In the cosmic variance limited 
case, because of = 1, S/N becomes a simple increasing function: {S/N)bb = 

(f'max + 3)(£max ” l)/2 (black dotted lines), independently of the shape of Cf^ and the 
values of input model parameters. The S/N is lower in the Planck-\\ke realistic experiment; 
however it can exceed one in the examples with the smallest values of A:* shown. A greater 
BB signal can be obtained at larger values of 5. 

The sourced GW signal breaks parity, generating a nonvanishing correlation between 
the CMB temperature anisotropy and B-mode polarization (TB) [35-38, 57]. In Figure 7 we 
compute the forecasted signal-to-noise ratio for the detection of such signal: 



2 

TB 


^max 

^( 2 £ + l) 


e=2 


{cpf 

r'TTr'BB 

^£,dat 


(4.11) 


in the Planck-like realistic experiment and the ideal cosmic variance-limited experiment. In 
the Planck-like measurement, we can neglect the noise spectrum of temperature mode NJ'^, 
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Figure 8. First row: Forecasted signal-to-noise ratio for the detection of the sourced TTT bispectrum 
in a cosmic variance-limited, CMB experiment, as a function of the maximum multipole i included 
in the analysis. Second row: TTT reduced bispectrum for equilateral triangles of multipole 1. The 
colored lines are for the same parameters as the previous figures. Solid (dotted) lines refer to the 
contribution to TTT from the sourced scalar (tensor) modes. For comparison, the black dashed 
line indicates the forecasted S/N ratio for scale-invariant equilateral non-gaussianity of magnitude 
/rl"* = 150 (a rough indication of the level that can be reached by Planck in the scale-invariant case 
[58]). 


since this is negligibly small compared with the signal on our interesting scales {i < 500). 
We see that, among the examples shown, such a signature can be marginally detected only for 
the examples shown at (5 = 0.5 in the cosmic-variance-limited experiment. Not snrprisingly, 
these are the cases that also lead to a stronger BB signal, cf. Figure 6 . Among the examples 
we considered, these are the only cases for which Tpeak (namely the tensor-to-scalar ratio 
evaluated at the top of the GW signal) is substantially greater than 10“^. Based on the 
difference between the 5 = 0.2 and 5 = 0.5 cases, it is very likely that a large TB correlation 
can be produced at greater values of 5, possibly at a detectable level even in a Planck-\\ke 
experiment. As a comparison, we note that [37] forecasted that a one a detection is possible 
for r ^ 0.002 in the scale invariant case, and a few a detection for the largest values of r ~ .05 
that can be achieved in our scenario. The results of [37] are slightly more optimistic than 
ours because of the scale invariance of the correlator considered in that paper, that allows to 
use information from a larger number of multipoles. 

4.2.2 CMB bispectra 

Let us now study the possibility of detecting the non-gaussian statistics of the sourced modes. 
All the CMB temperature and polarization bispectra due to the scalar and tensor non- 
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gaussianities, discussed in the following analysis, are computed by means of the techniques 
outlined in [59, 60]. 

In Figure 8 we show the forecasted signal-to-noise ratio for the detection of the sourced 
TTT bispectrum as a function of the maximum multipole i included in the analysis, reading 


n) 


TTT 


E 


IrTTT I 


criTT r^TT r^TT 
£i,£2,4=2®^^! ^^2 


(4.12) 


with 




E 

m\ra2m^ 


h h h ^ 

mi m2 m3 y 


(0£l mi Oil2m2 ) 


(4.13) 


denoting the angle-averaged bispectrum. Here we show only the results in a cosmic variance- 
limited CMB experiment. The same results will be obtained also in the Planck-like experi¬ 
ment, since the instrumental noise is perfectly negligible for i < 500. All the lines shown are 
for the maximum amount of particle production allowed by the WMAP TT data (namely, 
for = ^*qimit shown in Figure 4). 

Given that the sourced signal has a much greater deviation from gaussianity that the 
vacuum one [41], and given the stringent limit on non-gaussianity, one may have expected 
to find detectable non-gaussianity in the examples shown. Figure 8 shows that this is not 
the case. The main reason is that, for the values of /c* that we have chosen, the sourced 
signal manifests itself only at the largest CMB scales, so that only those scales are relevant 
for the phenomenological study of the bispectrum. This is confirmed by the colored lines 
shown in the Figure, where the signal-to-noise ratio saturates at values of £max well below 
those included in the Planck studies (in particular, the smaller /c* is, the smaller is the value 
of f'max at which the signal-to-noise ratio saturates). This contrasts with the scale invariant 
case shown by the dashed line, where we can see that S/N continues to grow with increasing 
^max- The strong weakening of the non-gaussianity limits with decreasing values of ^max can 
be also seen for example in Figure 11 of [58]. 

In the analysis, we distinguish between the scalar tensor contribution 

to the TTT signal. In the tensor case, both even (£1 + ^2 + ^z = even) and odd (£1 +(.2+(-z = 
odd) are included in the analysis [57, 61]. We can also see all this from the second row of 
Figure 8, which presents the TTT reduced bispectrum, defined as 




(2G + 1)(2£2 + 1)(24 + 1) 


dvr 


h h h 

0 0 0 


n ~l 


B. 




(4.14) 


on equilateral triangles as a function of the multipole We observe that, contrary to the 
scale invariant case, the bispectrum significantly decreases at the higher multipoles shown 
in the figure (also in this case, the smaller fc* is, the sooner in t space the bispectrum starts 
decreasing). Also in the second row we distinguish between the contribution of the scalar 
and the tensor mode. While in the examples with 6 = 0.2 the scalar contribution dominates 
over the tensor one, the tensor mode plays a non-negligible or even dominant role in two of 
the examples shown at 6 = 0.5. 
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Figure 9. Left panel: Forecasted signal-to-noise ratio for the detection of the sourced BBB bispectrum 
in a realistic CMB experiment with Planck-\\kB sensitivity (Colored solid lines), and in a cosmic 
variance limited, CMB experiment (Colored dotted lines), as a function of the maximum multipole ^ 
included in the analysis. Right panel: BBB reduced bispectrum for equilateral triangles of multipole 
The examples shown are for 5 = 0.5; the other parameters are chosen as in the previous figures. 


Another interesting non-gaussian signature of the sourced tensor mode is the BBB 
correlation. The signal-to-noise ratio is computed as 



2 

BBB 


^max 

E 





c/^BB /^BB /~<BB 
,dat '-"£2 ,dat '-^£3 ,dat 


(4.15) 


Like the BB case, we here analyze the detectability both in a Planck-V\ke experiment and in a 
cosmic-variance-limited experiment. In Figure 9 we show that this signal has a much higher 
prospect of being detected than the TB one in the cosmic-variance-limited measurement. 
Such large S/N seen in the figure is due to the fact that no instrumental error is assumed, 
so that all the noise is due to the small BB signal (instead, the largest TT signal due to 
the vacuum scalar modes constitutes a stronger noise for the TB measurement in the cosmic 
variance-dominated case). 

We also computed the B-mode bispectrum for 5 = 0.2, obtaining a larger S/N ratio 
than in the 5 = 0.5 case shown in Figure 9. We do not show this result, as we are concerned 
about its reliability: the expression for the noise in (4.12) and (4.15) assumes that the signal 
is approximately gaussian. We investigate whether this is the case in Appendix F, where we 
estimate the departure from gaussianity of the full C and /i+ signals. Our estimate suggests 
that the departure is very small in the scalar case, and marginally small for /i+ at (5 = 0.5. 
However, this is not the case for 5 = 0.2. The result performed there suggests that the noise 
in this case may increase by a O (1) factor with respect to what we computed using (4.15). 
This still likely results in a potentially visible signal in a cosmic-variance-limited experiment, 
but, due to this uncertainty, we do not show this result here. 

As in the TTT case, the S/N shown in the Figure 9 is computed by summing all the 
signals in = even and li+l 2 + H = odd. The sum enhances the S/N. On the other 

hand, the analysis in the limited domain, ii + ^2 + t-z = odd in TTT or +12 + = even 

in BBB, is also very informative, since such signals originate from nonvanishing parity-odd 
non-gaussianity [62, 63], namely a distinctive signature of this model. More comprehensive 

^®See [57, 61, 63-66] for the other possibilities to generate parity-violating CMB bispectra. See also [58, 67, 
68 ] for the WMAP and Planck constraints. 












analysis with the other auto- and cross-bispectra including the E-mode polarization, not 
discussed above, may also help to improve the detectability [56-58, 69-71]. 

5 Conclusions 

Given the expected improvements in sensitivity of the measurements aiming at the detection 
of primordial tensor modes, it is crucial to determine whether there are viable alternatives to 
the standard mechanism of creation of GW through the amplification of vacuum fluctuations. 

In the present paper we have shown that, if an axion-like held that is only gravitation¬ 
ally coupled to the inhaton experiences a transient epoch of relatively fast roll, then tensor 
perturbations can receive an additional contribution whose amplitude is not proportional 
to the energy scale of inhation. This implies that, for a broad range of parameters, ten¬ 
sors of inhationary origin might be detectable even if inhation happens at a low scale, and 
even for subplanckian inhaton displacements, without contradicting the existing constraints. 
Remarkably, the simple potential (1.2) naturally satishes our requirements. To leave an ob¬ 
servable GW signal, the transient epoch must occur while the large scale GMB modes leave 
the horizon. 

The system we have considered comes with a rich set of signatures. First, the amplitude 
of the induced tensor modes can be strongly scale dependent. As a consequence signihcant 
R-mode GMB huctuations might be detected only at scales related to the reionization bump, 
or only at the B-peak at £ ~ 80, but not at both angular scales. In particular, detection at 
£ ~ 80 without a corresponding detection at ^ ~ 10 wonld be interpreted as a locally blue 
tilt (similarly to the situation considered in [45]) for the primordial tensor spectrum, that is 
usually considered as a signature that would falsify the inhationary scenario. 

Moreover, the parity-violating nature of the process generates a characteristic {TB) 
correlator, that, in an ideal experiment, might be marginally detectable if the amplitude of 
the induced tensors is large enough. 

Finally, the statistics of the produced gravitational waves has a departure from gaus- 
sianity that is signihcantly more marked than for the vacuum modes. We have shown some 
choices of parameters for which a cosmic variance limited experiment could detect a nonva¬ 
nishing (BBB) correlator at a very high level of signihcance. Moreover, we have computed 
the phenomenology of the model only for a very few choices of parameters. As we mentioned 
in the previous section, a larger tensor mode can be produced for increasing values of the 
axion mass (greater <5), and it is likely that this can lead to larger TB and BBB correlators, 
which may be observable (particnlarly in the BBB case) already in a P/ancA:-like experiment. 
A more comprehensive analysis with the other possible auto- and cross-bispectra including 
the E-mode polarization may further enhance a signal-to-noise ratio. We believe that our 
result motivates a dedicated study of such localized signal in the Planck data. 

For the examples we have studied the effect on the temperature correlators is statistically 
small. This is the case because (i) due to helicity conservation, the amplitude of the sourced 
scalar perturbations is smaller than that of the tensors, and (ii) all relevant effects take place 
at multipoles with £ < O(IO^), where cosmic variance is large. We have thus explicitly verified 
that the model studied here is one of the very few existing examples for which a sourced GW 
signal is produced, while respecting the strong limits imposed from the observed T modes of 
the GMB. 

On the other hand, different choice of parameters can be made in the same model that 
can resnlt in a stronger scalar mode production than a tensor mode. This can easily be the 
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case if the field a rolls for a considerable amount to e-folds during inflation. This would 
generate a bump in the scalar modes with significant departure from gaussianity and with 
a greater width than in the examples considered here. By tuning the moment at which the 
roll of a takes place, the bump can manifest itself either at CMB or LSS scales. 

Remarkably, all the above features originate from the single operator appearing as 
the last term of the lagrangian (2.1), that is naturally expected to occur in theories with 
pseudo-scalar fields. It would be interesting to perform a more comprehensive study of the 
scenarios involving just Einstein’s Relativity and a local quantum field theory that can lead 
to comparable phenomenology. 

The next few years will witness major improvements in the observations aiming at the 
detection of primordial B-modes. It is important that all potential signatures are considered 
when the new data will arrive. 
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A The expression of A+(r,/c) in WKB approximation 

In this appendix we derive the expression (2.II) for the mode function of the photon. We 
start from the equation of motion (2.10) for the positive helicity photon, that, after defining 
X = —kr, X* = reads 

d^A+ ^ A _ 2 2g, 

\ X + {x/x^)~^ 

and that we solve using the WKB approximation. To do so, we write it in the general form 


A+ = 0, 


(A.I) 


A'Ux) -I- Q{x) A+(x) = 0, 


(A.2) 


where Q{x) vanishes for x = x and is positive for x > x, so that we can define 

g(x) ^ ( YY- ’ Y - 

( —k{x) , X < X , 

Then, for x > x, in the adiabatic regime |y(x)| <C p{xY, the general solution reads 

^ (i ™ (i " i)' 


(A.3) 


(A.4) 


with a and f3 constants, and where the symbol = denotes approximate equality in the WKB 
approximation. Similarly, for x < x, the solution in the WKB regime |Y(x)| <C k{xY 
will be a linear combination of exp ( — K{x)dx j / k{x) and exp ( K{x)dx j /yY(x). To 
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determine the coefficients of such a combination in terms of a and /3 we expand eq. (A. 2 ) 
near x = x as 

A'^(x) + Q'(x) (x — x) A+(x) ~ 0, (A.5) 

that can be solved in terms of Airy functions, and that we join to the adiabatic solutions in 
the regimes x — x —?■ ±oo, as discussed for instance in chapter 7 of [72]. The final result is 


A+(x < x) = 


a/2 

\/k(x) 


exp I — / k{x) dx + 


/3 




exp 


k{x) dx 


(A.6) 


The requirement that the photons are in their adiabatic vacuum in the ultraviolet trans¬ 
lates into the boundary condition a = regime x <C x we keep only 

the term that corresponds to the growing mode, so that 


A+(x < x) ~ — 


\/2k k{x) 


exp 


n{x)dx 


(A.7) 


As a check of the validity of our procedure we can consider the limit <5 —)• 0, so that 
2 B 

x^+x-^ brings us back to the case of a constant velocity for the rolling axion hrst 

studied in [32], Indeed, in this regime 



1 dx 


2 ^ arccos(y^) - a/x ( 2 ^ - x) ~ - 2 a/^ 


X . 


(A. 8 ) 

(where in the last step we have expanded to first order in ^/x), that, once inserted into 
eq. (A.7), gives a result that coincides with the approximate expression given in [32]. 

In the general case (5 7 ^ 0 it is impossible to compute the integral of k{x) in closed form. 
Nevertheless, we can write 


px px px 

/ k{x) dx = / k{x) dx — / k{x) dx , 

Jx Jo Jo 


(A.9) 


so that only the second term depends on x. Since we care only about the regime of moderately 
small X, we consider only the leading part of the integrand as x —)• 0: 


px px 

/ k{x) dx ~ / 

Jo Jo 


4 ^* 


X (x/x^ 




dx = 


4^ 

1 + d 


1+1 
X 2 



(A.IO) 


This equation, once inserted into eq. (A.7), gives the full dependence of the function A_|_(r, k) 
on the conformal time r as it appears in eq. (2.12) in the main text. The overall normalization 
of A+(r, k), on the other hand, has to be computed numerically, and we determine it by 
matching the WKB approximate solution with the numerical solution of eq. (2.10), yielding 
the constant iV(^*, x*, d) defined in eq. (2.12). We have checked that, in the part of parameter 
space we are interested in, the expression found this way approximates the exact solution at 
the 30% level or better. 


B Sourced scalar modes at constant 

In this Appendix we show how, for constant values of the isocurvature perturbation as¬ 
sociated to a is partially converted into curvature perturbation during superhorizon evolu¬ 
tion [47, 48]. As a consequence, large fluctuations in the gauge mode can introduce a signif¬ 
icant sourced component in the curvature perturbations. We start from the equations (3.6) 
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already given in the main text, and we redefine 


)Vf2 = _ 2y^>2 ^ + 9e0 + Sec. - 3?7^ 6,/e^ 

2 + 9e, + 3e^-3r?, 

If Mf- is approximately constant we can diagonalize the mass matrix as 


(B.l) 


M = C^AC, C = 


cos 0 sin 6 


— sin 9 cos 9 

so that we can solve the system (3.6), obtaining 


A = ( -^0 0 

' 0 A. 


(B.2) 


= sin 9 cos 9 / dr' 


G^nr,r')-G^(r, r') sJT',k) , 


(B.3) 


where (r, r') is the retarded Green function associated to the operator + ^ 


Gki'T, t) = 0(r - r) {J^^i-kr) Y^{-kT) - Jf^i-kr') Y^{-kT )) , /i = \j^ + \ ■ 

(B.4) 

Inspection of the equations given in subsection 2.1 above shows that, for ^ > 0(1)) 
the source term is exponentially suppressed for |fer'| > 1. As a consequence, we 

can set |A:r'| <C 1 in the propagator. Also, since we are looking at modes far outside of the 
horizon, we have |A;t| <C 1, and we can use the approximate expression for the propagator 


Ofc(r, t ') ~ 




(B.5) 


To first order in the slow roll parameters, denoting = 2 + 5\a and = 2 + we 
can expand 

' 


G^(r, r') 


'TT 


1 - ^ ^ log — 

3 3 T 


(B.6) 


so that, neglecting term —(IA/3 which is subdominant with respect to ^ log 


~ sin^ cos^^l^^-y^ J dr' ^log GI{t, T')Sa ■ (B.7) 


Since most of the integral in dr' gets contribution from t ' ~ — 1/fc, we can approximate 
log ^ ~ — log(—/cr) ~ Nk, the number of efoldings of inflation from the time the mode k 
has left the horizon. 

As a consequence the perturbation in 8(j) is the same one that would have been obtained 
in the case in which a = (j) (first studied in [41]) times the factor sin0 cos 9 _ 

‘2y/GrGj> Nk, that accounts for the fact that in this model a ^ (f>, and that quanta of cj) are 
obtained from the conversion of quanta of a. Our computation reproduces the result of [47], 
up to a factor of 2. The origin of such a difference can be traced to the second equality 
in eq. (5.17) of [47], where a is treated as a constant. By properly accounting for the time 
dependence of d = y/Te^ H Mp, the quantity 2e^ —A 2 in the last term of that equation should 
read 3e^ — X 2 , leading to a result equivalent to ours. 
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C Scalar Modes 


In this appendix we show the explicit derivation of the scalar power spectrum and bispectrum 
in the model ( 2 . 1 ), whose fitting functions used for the phenomenological studies are given 
in (4.7). 

Starting from (3.11), we split into the vacuum mode and sourced contribution in the 
way as in (3.12). Using the solution (3.15) for we find the sourced mode of C 

^ j dr" Gkir',T")Sa{T",k) . (C.l) 

where and Gk{T,T') are defined in (3.10) and (3.16), respectively. The explicit expression 
for Sfj is obtained by using (2.15) as 


same 

as 


(t, 


Sa{T,k) = a— 


3 , ^3 




f J (27r)3/2 

r d^p (+) j+) 

4/ J (27r)3/2 


(p) e-"'" (fe — p) p^^^\k — + \k — 


xA{T,p)A(T,\k-^) d+(p) + d^: {-p) d+{k-p) +d'!:(-k + p 


A 


(C.2) 


after symmetrizing p and | /c — p|, where 

A{t,p) A(T,\k - p^ = N[i^,-pn,6]N i*,-\k - p\t^,5 


X exp 


4^y2 . X./2 


\f-pT+G-\k-p\T 


1 + (5 \T, 

Since we are interested in large-scale modes, we can safely assume —kr <C 1, leading to 


(C.3) 


Gk (r, r') ~ 0 (r - r') '^ 3/2 {-kr') , -/cr < 1 


(C.4) 


while this approximation is not valid for Gk plugging this and (3.16) into (C.l), we 

have 


^r, itj ~ 


37r3/2F r dr' 


2M„ A;3/2 /__ r 


— J'i/2{-kT') Vecr(r') / dr'W- t" Sair" ,k) 


h /2 [-kr') y3/2 {-kr") - ^3/2 {-kr') J 3/2 {-kr") . (C.5) 


We then plug (C.2) into the above expression and rescale 

P 


p = - , x' = —kr' , x" = —kr" , 
k 


(C. 6 ) 


to obtain 




8M„/ 


dAp 

(27r)3/2 ■ 


+ \k — 


xN , p X*, (5] IV \k - p\x^, 6 


1 


■ 

V 

p, k 

Tc 


C*, X*, d, y^+\J\k-p\ 


(C.7) 
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where we have defined 
p,k 


V 

Tc 




«+(^ + a+(-^ a+[k-^+a\[-k + ^ 

4^y2 ^//(l+5)/2 


- r. e _ /■“ dx' , e^ix') , 

c[^*i x^, 5, Q] ^ — J 3/2 (x ] W- / dx 

Jo X W €(7^* J 


n 2;"3/2 gjjp 


l + <5 


5/2 




X [-^ 3/2 ( a^O ^ 3/2 ( a ;") - ^ 3/2 ( a ;') J0/2 ( a ^")] > 


(C. 8 ) 


where we have sent —kr —)> 0 for the lower bound of the integral. The functional forms of ^ 


and €„ is shown in (2.6). Using this, and the fact that 

2?* 


^CT - 


i{x) = 


y/2MpH 

2 


= const. X we have 


{x/x,y + (x*/x)' 


l egjx) _ _ 

Co-,* (x/x*)*^ + (x*/x)^ 


(C.9) 


where x* = —/cr* = [k/a{T^)] /H is the ratio between the physical momentum of the mode 
and the horizon at the moment in which ^ is at its maximum. We emphasize that apart from 
the creation and annihilation operators, the dependence on k arises only through x* = k/k^ 
in the expression (C.7). Using (C.7), we proceed to the calculations of the scalar power 
spectrum and bispectrum in the following subsections. 

C.l Scalar Power Spectrum 

The power spectrum of curvature perturbations C, can be defined as 

Pc {k) {k + k')^^ (C (k) C (fc') ) , (C.IO) 

and it consists of vacuum mode and sourced contribution, namely. 


Vc{k) = vy’{k)+vy\k). 


(C.ll) 

m. 


Notice that these two modes are uncorrelated, and thus there is no cross term. Here (fc) 
is the standard vacuum mode, for which, using eq. (3.14), we find 


V^^\k) =_ — _ 

^ S'rr'^e^M^ ’ 

neglecting the slow-roll corrections. For the sourced mode, we compute, using (C.7), 


(C.12) 




dTT^H^a^e, 


CT,* 


d^p d^p' 


UM^P J (27r) 

1/4 


V 


p,k V P,k' 


X (j)p\k — p\\k' — p\j + \k — + \k' — 

x/V px*, (5] Ai P,\k—p\x^,6 N [P, p'x^,, 6] N |U — p'| x*, <5 


xTc 


C*, d,^/p+y/\k^\ Tq C*, d, P, y/f +\J\k' 

A'* A'* 


, (C.13) 


where k^ is the mode that exits the horizon at the moment when ^ takes its maximum value 
namely fe* = a(r*)77. Evaluating 


V 


p,k 


V 


p', k' ) = \^k + kj {p + p') + \^k — p + p'^ 

{p) {k-pj (y') (fe' -p') , (C.14) 
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and using the definition (C.IO), we obtain 




d^p 


U ^ 64M2/2 J (27r)3 


4+' (pI 4*’ 


[C*, px^, 6]N‘^ \k-p\x^:,5 


k — p] p^^^\k — + \k — 


T/ 


C*, x^, 5,y/p+ \J\k-p\ 


(C.15) 


Expressing H^/Mp and a^eo-,*//^ in terms of and respectively, by (C.12) and 
= &'^{n)/{2M^H'^f‘^) = we find 


vf{k)^ /2,c(e*,x*,5) 


,(o) 


1 2 


(C.16) 


where 


/2,(^(C*5 3^*? o) — O ^* / /'O AS I ■" ~ it' ~i 

2 J (27r)3 \ p |A;-p| 


^ p k — p 


+ \k — 


xN"^ px^, S\N‘^ i^,\k-p\x^,5 


T/ 


c*, X*, 6,y/p+ \/\k-p\ 


(C.17) 


using the property of the polarization vector, (p) ^ = (1 — AA'p • /4. For the 

concrete evaluation of the p integral, we denote the cosine of the angle between p and k by 
r]. After taking the trivial angular integral, we have 


/2,c a:*, (5) f dp [ dTjp^^"^ {l-2pr] + p‘^)^^'^\p^/^ + {l-2pr]+p‘^y^^^ 


1-1 


1 2 


1 + 


p-rj 


(1 — 2pr] + 


[^*, px^, SjN"^ C*, (l - 2pr/+p^)^'^^x*, 6 


xTr 


C 


, X*, (5, p^/^ + (l - 2p?7 + p^) 


■ 2 ^ 1/4' 


(C.18) 


Alternatively, we can change the variables of integration from p and p to x and y such that 
X = p+\k — p\ and y = p — \k — p\. Then the integral reduces to 


Qvr^ „ r°° f 


(1 - yx + y + y/x-yy 

dy - 


Vx + yy/x - y 


X 


. x + y 

s*) 2 Xjf, 0 




. x-y 
s*) 2 x^, 0 






Vx + y + Vx-y 


V2 


(C.19) 


We evaluate /2 numerically using either (C.18) or (C.19). The total power spectrum of the 
curvature perturbations is, from (C.ll) and (C.16), 




(C.20) 


where P® is given in (C.12). 
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C.2 Scalar Bispectrum 

We define the bispectrum of the curvature perturbations as 

(ku k 2 , h) 5(3) (ki + k 2 + h) = (c (h) C (^ 2 ) C (^ 3 )) , (C. 21 ) 

where 13^ depends only on the magnitudes of the three momenta, under the restriction that 
they form a triangle. The bispectrum in principle consists of the vacuum and sourced modes, 
similarly to the power spectrum. However the vacuum mode bispectrum is small, and 
we focus on the contribution from the source, i.e. 

B^ ^ . (C.22) 


Taking the 3-point function of in (C.7), we have 

337r9/2f/'6Q,3g3/2 n ^3^^ 


3 

i=l 

xTc 


V 


Pi, h 


V 


P2,k2 


V 


29M3/3 J (27r)9/2 

+ \ki piX^, 5]iV \ki -Pi\x^, 6 


P3,k3 


3/"^ /- 1/2 


C*, Y' + \/l^* - 


Pi\ 


(C.23) 


Evaluating the vacuum expectation value, we find 

f d^p 


T. E\ 3V/^-H‘aY 1 

B}. ( / ci , k2, / C 3 ' 


26M3/3 kfk^k^J (27r)9/2 


p 


p,p + ki,p- ks 


p\p + ki\\p- ksl 


X (^Vp + \/\p + h\^ (\/\p + ki\ + \/\p-(\/\p- M + Vp 


xN^ 


xTc 


xTc 


^ P jr 




^ \P + ki\ , 

7 ? ^ 


/c* 




C \p-k3\ ^ 

S*? 7 ? ^ 


fc* 


^ ki Vp+y\p + ki\ 

TIT 




c h 

s*? 7 ? 


P + + \ P-h 


Vh 


c h s 
?*) , ) 0, 
'bjk 


p- A:3| + Vp 

Vh 


(C.24) 


where 

P [vi,V 2 , V 3 ] = (Ti) eS+C (i;2) (i;2) (i;3) (T3) (vi) 

VI-V2+V2-V3 + V3-V1 + {vi ' hf + (1)2 ' 1 ) 3 )^ + {h ' Dl)" 

O - 

+ (Di • h) {v 2 ■ h) + {h ■ V3) in ■ vi) + (D3 • Di) (Di • h) - (Di • h) (D2 • h) (D3 • Di) 


+ 0 vi ■ ih X V3) {1 + vi ■ V2 + V2 ■ h + h ■ Vl) 

O 


(C.25) 
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We can disregard the imaginary part as the bispectrum is real (see Appendix E). Rescaling 

u-i -^2 _ fes _k p 

fe = fci, X2 = -r , X3^ — , = — , P = Y , (C.26) 

k k k^ k 

and using the relations 


we obtain 


m! 




( 0 ) _ \/ 2 ^* 


/ Mp ’ 


( 0 ) 


1 3 


where 

/3,c (^*,a;*,(5,X2,X3) = 


T2T2T2 /3,C 3^25X3) 

/C2'23 

d?p 


(C.27) 


(C.28) 




2^2 


XoX 


2-^3 


( 27 r) 9/2 


(' 


^{p 


p,p + ki,p- x^k^ 


^\lp\p + ki\\p - X3k3\ y\/p + \l\p + [\/\P + ^A + \/\p-X3k3\j [\/\p-X3k3\ + y/p 

xN'^%, px^, 5 ]N'^ |p +fell X*, (5 iV^ I*, Ip - a^sfesl a;*, ^ 


X*, ( 5 , v^ + V Ip + fell 


xTr 




P + fell + Y |p - X3k3\ 

\7^ 


Tt 


C* 7 a^sa^* 7 ^7 


|p - a;3fc3| + y/p 


(C.29) 


where 3? denotes the real part. 

In order to evaluate the integrals, we can orient ki along the x axis and express /c 2 and 
k 3 in terms of X2 and X3, namely, 

ki = k{l, 0 , 0 ) , 

fe2 = ^ (-1 - a;2 + xl, y/- (1 - X2 + X3) (1 + X2 - X3) (1 - X2 - X3) (1 + X2 + X3), 0 ^ , 

fea = ^ (^-1 + 3^2 - xl, -\/- (1 - X2 + X3) (1 + X2 - X3) (1 - X2 - X3) (1 + X2 + X3), 0 ^ , 

(C.30) 

and then perform numerical integration. The phenomenology of this result is extensively 
studied in Section 4. 


D Tensor Modes 

In this appendix, we present the detailed derivation of the tensor power spectrum and bis¬ 
pectrum from the source, parallel to the scalar counterparts derived in Appendix C. The 
phenomenological study in Section 4 is done with the use of fitting functions (4.7) for the 
results derived in this appendix. 

We define the tensor mode operators in the basis of circular polarization, 

feA (r, fe) = A hij (^T, tj = ^ Qx (^r, tj , (D.l) 
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where the tensor canonical mode Q\ and the polarization tensor are defined in (3.19) and 
(3.20), respectively. The equation of motion for the canonical operator Q\ is shown in (3.21). 
We can decompose the solution into the homogeneous and particular ones, corresponding to 
the vacuum and sourced modes, respectively, in the same way as for the case of the scalar 
perturbations (3.12). The solution for the vacuum mode is given in (3.22), and the sourced 
contribution can be solved formally by (3.23). This, together with (D.l), provides the 
solution for the sourced mode as 

Gk{T,T')Sx{T', (D.2) 

where we have used a(r) ~ the Green function is given in (3.16), and the source 

5a is defined in (3.21). The explicit expression for 5 a can be found by using (2.14) as 


5a [T,k) = -—Tlij,x[k 




\/2Mp 


(27r)^/^ 

dx’p 

(27r)^/2 


® ^-ik-x 


Vx 


Ei Ej + Bi Bj 


k,p,k — p 


p 


1 + 


—r 


2 ^(t) 




xA{t,p) A[T,\k-p\] dx{p) + d\{-p) dx[k - p] + d\[-k + p 


(D.3) 


where and dx are the creation and annihilation operators, respectively, of the helicity 2A 
graviton modes. Recalling the relation (3.20), we have defined 


Vx 


k,p,k-p 


= 


k ) (p) ( k ) [k - p] . 


d+) 


Since we are interested in large-scale modes —/cr <C 1, we can approximate 

Gk (t, r^) ~ ^ [fer' cos {kr'^ — sin {kr')] , —kr <C 1 , 

which is identical to (C.4). Combining (D.2) with (D.3) and (D.5), we obtain 


^x'' ( D ^ ^ 


v/2F2 f d^p 


M^kV^ J (27r)3/2 


xN 


Vx 
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t X 
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u \k-p\ 

, u, ^ ^ 


ttA (p) + 4 (-^ 


dx[k-p) +d[[-k+p 


where 


Th[^*, x^, 6, ,p, q] =Thi x*, <5, v^+ a/^ 




(D.4) 


(D.5) 


(D. 6 ) 


Thi x^, 6, Q] = j dx' {x' cos x' — sinx') ^ 


— exp 


i*, x^, d, x/p+V^ 

4^y^ 3,/(i-r5)/2 


Th 2 X*, 5, Q] = / dx' [x' cos x' - sinx') W 77 ^ exp 

Jo V ) 


1 + 6 

4ey'x'(i+^)/2 
'1 + <J xf/^ 
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■ (D-7) 
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Here we have rescaled the variables x' = —kr', and x* = — fcr* = [/c/a(r*)] /H is the ratio 
between the physical momentum of the mode and the horizon at the moment in which ^ is 
at its maximum. For the evaluation of the integral, recall 


^(x) 


2 ^. 

(x/x*)”^ + (x*/x)^ 


(D.8) 


Using the expression (D.6), we compute the tensor power spectrum and bispectrum in the 
following subsections. 


D.l Tensor Power Spectrum 

We define the power spectrum of the tensor mode as 

Vx{k) {k + k') = ^ (hx (fc) h' [k') ) , (D.9) 

and separate it into the vacuum mode and the component sourced by the gauge field, as 
in (4.2). The vacuum mode is the standard one, and we focus on the computation of the 
sourced contribution in this appendix. Using (D.6) and taking the two-point correlator of 
, we find 


(fe) {k'y^ ~ {k + k'^ k,p,k-p Vy k,p,k-p \/p\k-p\ 


xN^ 


c P X 






q-2 

I h 


p \k-p\ 
s*? ^ k 


One can show that 


dc^Vx 


k,p,q 


VI, 


k,p,q 


= ^/#(l + Afc-p) (^1-t-Afc-g^ , 


. (D.IO) 


(D.ll) 


and therefore we find the sourced part of the tensor power spectrum 


pW (k) cx 




f d?p 

k^ J {2 tt)^ 


l + Xk-p]^ ll + Xk- ^ ^ 


\k-p\ 


\/p\k-p\ 


xN^ 


f P X 




^ X 

^5k 


fh 


p \k-p\ 

S* ? ; D, ^ ^ 


After rescaling p = p/k, we arrive at the expression 


pW (fc) f2,x{C.,x,,5) , 

where V^^^ is given in (C.12) and 

2 / d^p 


,(o) 


1 2 


p(o) 


/2,A (?*, X*, <5) = Svr 


3 ^l +Afe-p)^ [l + Xk ^ ^ 


\/p\k-p\ 


(D.12) 


(D.13) 


(2vr)3 \k-p\ 

xN‘^ px^, 6]N‘^ |fe-p| X*, 5 Th x^, 6, p, \k - p\ (D.14) 
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For the concrete evaluation of the integrals, we denote by rj the cosine of the angle between 
p and k. After taking the trivial angular integral, we have 


poo rl p5/2 (^1 _|_ _|_ \^l _ 2pr] + p^') 

f 2 ,x (?*, X*, 5 ) = 2 / dp dp - 


0 J-i (1 — 2jir7 + p2)3/4 

px*, (5] , x/l- 2pp + X*, 5 Th C*> ^7*, <5, p, x/l- 2pp + p^ . 

(D.15) 


Alternatively, as is done for (C.19), we can change the variables of integration from p and p 
to X and y such that x = p + |A: — p| and y = p — \k — p\. Then / 2 ,a takes the form 


/2,A (C*, X*, 5) = I ^ dx 


(l — y"^)'^ (1 + Ax)^ 5 
dy^ — , ^ \ -N^ 


y/x + yy/x - y 


f ^x 5 

Ss* ? r* 5 ^ 


xN^ 


c x-y 

S*) rv X*, 0 


q-Q, 

Ih 


1 ^*, X*, 6 , 


x + y X - y 


(D.16) 


We can now readily evaluate / 2 ,a numerically. All the phenomenological features of the tensor 
power spectrum in this model are captured in / 2 ^a- Since we are choosing ^ > 0 (i.e. cr > 0) 
so that the A+ modes are produced through the tachyonic instability (see discussion around 
(2.11)), only the + helicity state of the tensor modes is efficiently sourced by the gauge field, 
leading to the hierarchy / 2 ,+ ;:?> / 2 ,-. 

The total tensor power spectrum for each helicity state is 


rx = 




7r2M2 


1 + 




Gdvr^M^ 


/2,A (?*, X*, 5) 


= 8 esV, 


( 0 ) 


1 +v^f^/2,A (e 


8 


<5) , (D.17) 


and the tensor-to-scalar ratio is, using (C.20), 


r = 


Ea^a 


l + f|^f^(/ 2 ,+ + / 2 ,-) 
l + 62 pW/ 2 ,C 


16e. 


(D.18) 


Interesting phenomenological features arise when the sourced contribution dominates the 
tensor spectrum while it is subdominant for the scalar spectrum. 


D.2 Tensor Bispectrum 

We define the tensor bispectrum as 

^AiAaAg (fci,/c2,fea) (5^^) (^ki + k2 + fca) = (hxi hx^ {ko^ hx^ ’ (D-19) 

where in general BX 1 X 2 X 3 is real and depends only on the magnitude of the three momenta, 
forming a triangle. While it consists of the two contributions from the vacuum fluctuations 
and the source effects, the former is unobservable, and we focus on the latter. Moreover, 
since only one polarization state of the gauge field is enhanced, the produced gauge quanta 
efficiently source only one of the tensor helicity states. We thus define and have 

4'^=^iAA = ^AAA, (D.20) 
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where the superscript (1) denotes the sourced mode. We compute the 3-point correlation 
function of with the same A using (D. 6 ), 


^ ^ (^ 2 ) ( fca ) ) ~ <5^^^ {ki + k2 + k3 
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(D. 21 ) 


(D.22) 


and recalling (C.12), we obtain 
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A^^ - \ 2 u 2 A hx ih, X,, 6, X 2 , X 3 ) 
^ 1 ^ 2 ^3 
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(D.23) 


where 


227/2^6 ^ 
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In order to perform explicit evaluations of the integrals in (D.24), we fix ki along the x axis 
and write /c 2 and fcs in terms of X 2 and X 3 , namely, 

ki = k (1, 0, 0) , 

fc2 = - (^-1 - + -s/- {I- X2 + X 3 ) (1 + X 2 - X 3 ) (1 - X 2 - X 3 ) (1 + X 2 + X 3 ), 0^ , 

^3 = ^ (-1 + xl- -a/- (1 - X2 + X3) (1 + X2 - X3) (1 - X2 - Xz) (1 + X2 + X3), 0^ , 

(D.26) 

and the polarization vector for a given momentum q can be written in terms of its components 
as 


d±) (r-A = 




qxqzTiqy^ql + ql + ql qyqz^iqx^Jql + ql + ql ^Jql + q, 


^ V J^x + qlJql + ql + q'^z Jql + qlJql + ql + ql Jql + ql + ql 


(D.27) 

Using these explicit forms, we can evaluate (D.24) numerically for any given set of parameters. 
In the present case with ^ > 0 (i.e. ct > 0), the positive helicity state of the gauge field is 
produced, and consequently only the positive tensor mode is efficiently sourced, resulting 
in / 3 ^_|_ » / 3 ,-- The phenomenology of the tensor bispectrum is featured by / 3 ^a and is 
discussed in detail in Section 4. 


E Shape and properties of the bispectrum 

Let us first verify that our expressions for the scalar and tensor bispectrum given above are 
symmetric under the exchange of the external momenta ki, and real. In this discussion, 
we refer to the scalar bispectrum for definiteness. However, all our statements also apply 
to the tensor bispectrum. To show that the bispectrum is invariant under the exchange of 
any two momenta it is enough to show that {ks, ki, k 2 ) = {ki, k 2 , ks) and that 

B^^'^ (fci, ks, k 2 ) = {ki, k 2 , k^). To verify the first equality, we express B^^^ ki, k 2 ^ 

through expression (C.24) (namely, we replace /ci —)■ fcs, /c 2 —)■ ki and k^ —)■ k 2 both at the 
left and right hand sides of eq. (C.24)), and we relabel the integration variable as p —>■ p — k^. 
Using the fact that the three external momenta add up to zero, it is immediate to verify that 
the resulting expression coincides with the expression (C.24) for B^^'^ {ki, /c 2 , k^). To verify 

the second equality, we again start from the expression (C.24) to express B^^'^ (jti, k^, k 2 ^ and 
we relabelp —)• —p—ki. Also in this case, we recover the expression (C.24) for {ki, k 2 , k^). 

To verify that B;^ is real we instead use the reality of C (x), and the consequent identity 
C = C* (“^)’ write B(^ (jti, k 2 , ks'j = ^i, —k 2 , —k^^. We then perform a 180° 

rotation around the axis perpendicular to the plane defined by the three vectors ki (recall 
that they define a plane since they add up to zero). Under this rotation, the three vectors 
change sign. Due to statistical isotropy, the rotation does not change the bispectrum, and 
so B'^ (^—ki, —k 2 , —ks^ = B*(- ^fci, k 2 , k^. From the two identities that we have just written 
we see that the bispectrum is real. 
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As mentioned, identical properties apply also for B\, due to statistical isotropy, and to 
the fact that h\ (x) is scalar and real. 

We just proved exact identities, that apply to our explicit final results for and 
/s.A- We now provide an approximate expression for the bispectra, written as a sum of 
factorized terms. Namely, each of these terms appears as a product of three functions, each 
of which depends on one and only one of the external momenta. As well known, a factorized 
approximate expression for the bispectrum considerably speeds up its use in the data analysis. 
We conjecture this expression by noting that we expect the bispectrum to have a bump at 
the value of momenta ki k 2 — — k^, at which the particle production is maximum, 

namely at which also the two point function is maximum. This leads us to expect a shape 
dependence of the type 

fsj {h, k 2 , ks) ~ Normalization x [f 2 j (ki) f 2 ,j (/C 2 ) f 2 ,j , (E-1) 


both in the scalar (j = Q and tensor (j = A) case. The normalization factor can be obtained 
by evaluating the bispectrum in the exact equilateral limit. A form that is factorized, and 
symmetric in the three momenta is 


h,j {h, k2, ks) ~ 


h,j (fci, ki, ki) hj {k2, k2, k2) h,j {ks, k^, k^) 

3^^ (fci) (fc2) 3^^ (fca) 

X [f2,j {ki) f2,j (/C 2 ) h,j , 


(E.2) 


The approximate expression (E.2) has the additional computational advantage that we simply 
need to evaluate /s in the equilateral limit. 

To verify the accuracy of (E.2), in Eigure 10 we show the comparison between the 
exact and the approximate bispectrum on isosceles triangles, k 2 = k^. The black/solid 
contours refer to the left hand side of (E.2), while the purple/dashed contours refer to the 
left hand side. The contour lines shown in the the figures are obtained by evaluating the 
bispectra on a grid of values in the region shown, and by interpolation (for the functions 
/ 2 ,f and / 2 ,A used to produce the dashed lines we instead employ the fitting functions given 
in Subsection 4.2). In each panel, a black (purple) dot indicates the triangle for which the 
exact (approximate) bispectrum is maximum. In both cases, the two dots actually appear as 
nearly superimposed or as superimposed to each other, as the locations of the two maxima 
are nearly coincident. We see that indeed the bispectrum is maximum on an equilateral 
triangle of scale /ci = A :2 = /cs ~ few x A;* approximately equal to the scale at which the 
sourced power spectrum is also maximum. Most importantly, as shown by the figure, the 
approximate bispectrum provides a very accurate description of the exact one (particularly, 
close to the maximum). 


F Deviation from gaussianity 


In this Appendix we estimate the departure from gaussianity of the statistics of the sourced 
and h!'^ modes. Eor brevity we denote either of the two helds as /, and we distinguish 
between the two fields only at the end of the appendix, when we make explicit evaluations. 
We assume that / is approximately gaussian (to be verified by this computation) and, to 
estimate the departure from gaussianity we start from the conventional definition of the local 
/nl^^ parameter 


/ (^) — fg (®) + /nl 




(FT) 
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fc* 

Figure 10. Exact (black solid lines) vs. approximate (purple dashed lines; eq. (E.2)) expression for 
the scalar (left panel) and tensor (right panel) bispectrum, for isosceles triangles, /c 2 = fca, and for 
= 5 and S = 0.5. The black (purple) dot indicates the triangle for which the exact (approximate) 
bispectrum is maximum. The contour lines indicate the triangles for which the bispectrum evaluates 
to (from inner to outer, respectively) a fraction of 0.9, 0.7,0.5, 0.3,0.2 than the maximum value of the 
exact bispectrum. The red line at X 2 IX 1 = 1, indicates the equilateral triangles, while the lower slope 
X 2 = xi{2 line indicates folded triangles (smaller ratios are not possible for isosceles triangles). 


where fg is gaussian. The statistics of / is nearly gaussian if the first term in this decompo¬ 
sition dominates, namely if 

/n“7j(J)<1 « K^(A“£*‘)"(/(*)/(x))<1. (F.2) 

When this condition is satisfied, we are also estimating that the presence of a nonvanishing 
three-point function does not change the two-point function significantly (namely, we are in 
this way estimating the loop diagram where the nonvanishing three-point function enters as 
a vertex). 

In Fourier space, the decomposition (F.l) reads 

f (^) = fg (^) + /nT^ J ^^^^3/2 fg fg - p) ■ (F- 3 ) 

Under the assumption of small departure from gaussianity, namely that the first term domi¬ 
nates in (F.3), we obtain the three point correlation function 


B{h,k2,h) = 2V27r^/y^^l-^ 


V{ki)V{k2) 

U3 l3 
^2 


iP(fcl) V{k3) V{k2) ViksY 

^3 ^2 ^3 J 


(F.4) 


where the two- and three-point correlation functions are expressed, respectively, in terms 
of V and B as in eq. (4.1). The constant corresponds to a specific shape of the 

bispectrum, denoted as local non-gaussianity (the standard local non-gaussianity parameter 
is obtained from using ( in (F.3), and from rescaling —>• | ^^e present case, 

we promote to an effective momentum-dependence nonlinear parameter /nl (ki), by 
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Figure 11. Left (right) panel: parameter (fc) (r+ (fc)) for the greatest value of used in the 
previous figures, for S = 0.2, 0.5 and for = 10“®. 


inserting in (F.4) the power spectra and bispectra of the sourced functions and We 
evaluate this expression in the equilateral limit, where the bispectrum of our sourced signals 
is maximum (this likely overestimates the departure of gaussianity of the sourced signals). 
In this way we obtain 


/NL,i (fc) 


{k, k, k) 
(k) 


i = C, A. 


(F.5) 


To estimate the departure from gaussianity we write the analogous of (F.2) in momen¬ 
tum space, accounting for the momentum dependence of /NL,i- We thus obtain the condition 


„ r *,2 rdkk'^si{k,kk) _ 

Using the parametrization (4.6), we obtain 


dk ,,, ^ 

— ri{k)<l. 


(F.6) 


,,, k^^Bf{k,kk) 
n {k) = . _ I _o ' = 


72 7r^Vf{k) 


6 '15(0)3 r 2 


fO 'p 


■^3,C 

72n^{l+elvf ^ hx) 

,3 'p(0)3 r2 
^ 4 > ' C ‘'3,A 


_ f2x) 


3 , * = C , 

T , i = A . 


(F.7) 


In Figure II we show the quantity (k) (left panel) and r+ (k) (right panel) for the 
greatest values of considered in the previous figures, and for which the departure from 
gaussianity is therefore greatest. The quantity (k) is maximized at the bump of the sourced 
signal, k ~ few x fe*. We note that (fc) <C 1, while r_|_ {k) < 1. The integrals of the curve 
shown evaluate to 


R(^ {5 = 0.2, = 4.5) ~ 0.00098 , Rc{6 = 0.5, = 5.3) ~ 0.015 , 

R+ [5 = 0.2, = 4.5) ~ 1.6 , R+{5 = 0.5, = 5.3) ~ 0.35 . (F.8) 

The C signal is approximately gaussian, as the two-point function is dominated by the 
vacuum mode (which thus dilutes the non-gaussianity of the sourced mode). The tensor signal 
is mildly gaussian (we recall that the expression (F.5) likely overestimates the departure from 
gaussianity). This explains why in the examples we have studied in the main text the S/N 
ratio for detecting BBB is greater than that for TTT. It is possible that the large BBB signal 
can impact (at one loop) the BB correlator, thus inducing a greater two-point function (and 
that, in general, higher order correlators may increase the lower order ones). However, the 
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fact that i?+ ~ 1 suggests that this cannot change the order of magnitude of the result we 

have computed (nor the conclusion that the model can produce a visible BB signal, while 

being compatible with the limits from XT). 
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